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20.  ABSTRACT  (Continued) 

striated,  high-altitude  nuclear  plasmas,  has  been  difficult  to  characterize  in 
a  general  and  systematic  way.  The  Wideband  satellite  experiment,  with  its 
multiple,  coherent  frequencies  and  its  variety  of  observing  geometries,  has 
provided  a  wealth  of  data  upon  which  to  base  a  general  model. 

In  this  work,  careful  and  detailed  analysis  has  been  performed  of  VHF,  UHF, 
and  L-band  data  from  the  Wideband  receiving  stations  at  Ancon,  Kwajalein,  Stan¬ 
ford,  and  Poker  Flat  to  establish  the  statistical  behavior  of  scattered  signals 
over  a  wide  range  of  scattering  strength,  Fresnel-zone  size  (which  sets  the 
effective  distance),  and  scattering  geometry.  Histograms  of  intensity  and 
phase  were  compared,  by  means  of  chi-square  testing,  with  probability  density 
functions  derived  from  four  leading  signal-statistical  hypotheses:  log-normal; 
generalized  Gaussian,  of  which  Rice  statistics  are  a  special  case;  two- 
component,  which  invokes  the  log-normal  and  generalized  Gaussian  hypotheses  in 
different  regimes  of  the  signal  fluctuation  spectrum;  and  Nakagami-m,  which  is 
an  approximation  to  generalized  Gaussian  statistics  addressing  only  the  behav¬ 
ior  of  signal  intensity. 

The  analysis  disclosed  the  expected  efficacy  of  log-normal  statistics  in  the 
near  zone,  limited  to  conditions  of  weak  intensity  scintillation,  and  a  ten¬ 
dency  toward  generalized-Gaussian  statistics  approaching  the  far  zone,  limited 
to  a  wodulo-2-n  description  of  phase.  It  clearly  demonstrated,  however,  that 
log-normal  statistics  do  not  provide  a  generally  valid  description  of  inten¬ 
sity  scintillation  and  that  the  generalized-Gaussian  model  does  not  character¬ 
ize  the  overall  behavior  of  complex-signal  scintillation.  The  two-component 
model  was  found  to  be  generally  efficacious  for  describing  all  aspects  of 
signal  behavior  in  all  regimes,  providing  a  systematic  transition  between  the 
two  classical  models. 

A  strong  conclusion  has  been  reached  in  support  of  the  Nakagami-m  distribution 
as  a  general  descriptor  of  intensity  fading  and  of  the  normal  distribution  as 
an  indicator  of  phase  behavior,  the  first-order  statistics  of  which  usually 
are  altered  little  in  post-scattering  propagation.  At  present,  there  is  no 
underlying  signal-statistical  hypothesis  that  leads  to  both  joint  Nakagami- 
normal  behavior  and  an  accounting  for  correlation  between  intensity  and  phase 
fluctuations,  which  is  commonly  observed.  It  is  suggested  that  such  an  hypoth¬ 
esis  be  formulated  inductively  and  tested  empirically.  Meanwhile,  the 
Nakagami-m  distribution  can  be  used  confidently  to  describe  intensity  behavior 
and  the  normal  distribution  to  describe  phase  behavior.  For  a  full  description 
of  the  joint  statistics  of  intensity  and  phase,  at  present,  one  must  employ  the 
two-component  model,  which  requires  six  moments  for  its  specification.  A  joint 
Nakagami-normal  model  would  require  only  three  moments,  two  of  which  are  the 
commonly  employed  intensity  and  phase  scintillation  indices. 

The  analysis  reported  herein  also  elucidates  the  manner  in  which  radio-wave 
observations  and/or  communication  transmissions  of  finite  duration  impose  an 
effective  outer  scale  on  a  scattering  medium  observed  in  the  intermediate 
zone.  For  anisotropic  scatterers,  the  imposed  outer  scale  is  highly  geometry- 
dependent,  requiring  caution  in  extrapolating  experimental  results  to  disparate 
operational  situations.  It  is  demonstrated  that  presently  available  propaga¬ 
tion  theory  does  provide  the  tools  necessary  for  proper  extrapolation. 
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I .  INTRODUCTION 


A.  BACKGROUND 

The  statistical  behavior  of  an  electromagnetic  signal  that  has  propagated  through 
a  random  medium  has  been  investigated  by  workers  in  several  related  fields.  A  general 
solution  to  the  problem  is  of  interest  as  a  scientific  challenge  to  understand  wave- 
medium  interactions,  and  identification  of  specific  signal-statistical  characteristics 
has  been  sought  for  application  both  to  remote  sensing  of  various  propagation  media  and 
to  the  efficient  transmission  of  information  through  them. 

The  most  direct  theoretical  attacks  on  problems  of  signal  statistics  involve  moment 
calculations  for  the  distribution  of  one  or  more  observables,  such  as  the  electromagnetic 
wave's  quadrature  components  or  intensity  and  phase  (Hewish,  1951;  Bramley,  1955;  Bowhill, 
1961;  Yeh,  1962;  Uscinski,  1967;  Strohbehn,  1968;  Barabenenkov  et  al ,  1971;  Tatarskii, 
1971;  de  Wolf,  1972;  Fante,  1975;  Ishimaru,  1975;  Rino,  1976).  Applications  range  from 
radio  astronomy  (Salpeter,  1967;  Lovelace  et  al,  1970;  Valley  and  Knepp,  1976)  and  iono¬ 
spheric  physics  (Fremouw  and  Lansinger,  1968;  Crane,  1976)  to  atmospheric  remote  sensing 
(Beard,  Kreiss,  and  Tank,  1969;  Furuhama  and  Fukushima,  1973)  and  both  optical  and  radio 
transmission  (Ochs,  Bergman,  and  Snyder,  1969;  Gurvich  and  Tatarskii,  1975;  Fremouw  and 
Rino,  1978). 

Only  in  special  circumstances  (e.g.,  in  the  Fraunhoffer  diffraction  zone  of  a  medium 
containing  refractive-index  irregularities  with  an  effective  outer  scale)  have  sufficient 
moment  calculations  been  carried  out  (e.g.,  Mercier,  1962)  to  completely  define  the  first- 
order  distribution  of  the  complex  signal  (e.g. , ’Rician) .  In  other  circumstances,  only 
numerical  (Buckley,  1975)  and  empirical  (Knepp  and  Valley,  1978)  approaches  are  tractable. 
Moreover,  the  vast  majority  of  theoretical  work  and  much  of  the  experimental  effort  has 
centered  on  the  statistics  of  signal  intensity  (e.g.,  Rino,  Livingston,  and  Whitney,  1976) 
because  phase  calculations  are  difficult  and  phase  references  often  are  not  available  in 
experiments . 

In  the  work  reported  here,  a  direct  experimental  attack  has  been  made  on  identifying 
the  underlying  statistics  that  characterize  the  joint  behavior  of  the  intensity  and  phase 
of  a  (radio)  signal  that  has  undergone  narrow-angle  forward  scattering  in  a  remote  (iono¬ 
spheric)  region  containing  (real)  refractive-index  irregularities  with  an  effectively 
infinite  outer  scale  and  an  effectively  zero  inner  scale.  The  data  employed  were  obtained 
by  means  of  the  DNA-002  Wideband  coherent  satellite  beacon,  which  transmits  ten  mutually 
coherent  signals  ranging  in  frequency  from  VHF  through  S  band.  The  signals  were  recorded 
at  ground  stations  located  in  the  auroral,  mid-latitude,  and  equatorial  zones. 
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Historically,  there  have  been  two  heuristic  views  of  the  manner  in  which  signal 
perturbations  develop  due  to  interactions  of  a  wave  with  a  random  medium.  The  first  is 
based  on  the  mental  picture  of  a  wave  propagating  in  essentially  straight  paths  along 
which  the  phase,  4>,  is  randomly  increased  and  decreased,  differently  on  different  paths. 
Applying  the  central  limit  theorem  to  the  random  walk  undergone  by  the  phase  in  tra¬ 
versing  the  medium  leads  to  an  expectation  that  the  phase  will  be  a  Gaussian  variate  (e.g., 
Bramley,  1955).  A  natural  extension  of  this  construction,  which  is  particularly  con¬ 
venient  for  application  of  the  Rytov  (1937)  perturbation  technique,  is  to  suppose  that 
the  real  and  imaginary  parts  of  the  logarithm  of  the  signal  are  normal  variates,  so  that 
the  intensity  obeys  a  log-normal  distribution  (Strohbehn  and  Wang,  1972). 

Close  to  the  medium,  the  foregoing  "multiplicative"  (Strohbehn,  Wang,  and  Speck, 

1975)  model  (so-named  because  the  complex  phasor,  E,  representing  the  wavefield  is  per¬ 
turbed  by  multiplicative  effects  in  adjacent  slabs  of  the  medium)  is  satisfactory  even 
when  multiple  scatterings  lead  to  strong  phase  perturbations.  It  becomes  suspect  (Wang 
and  Strohbehn,  1974a),  however,  whenever  significant  fluctuations  in  wave  intensity,  I, 
develop  (even  for  moderate  phase  fluctuation)  through  the  effect  of  propagation,  because 
its  attendant  perturbation  technique  for  calculating  E  breaks  down. 

The  alternate  historical  view  leads  to  the  "additive"  (Strohbehn,  Wang,  and  Speck, 
1975)  model,  in  which  component  phasors  are  added  and  the  central  limit  theorem  is  applied 
to  the  real,  x,  and  imaginary,  y,  parts  of  E  rather  than  of  its  logarithm.  Far  from  the 
scattering  medium,  where  propagation  has  produced  a  uniform  distribution  (modulo  2 tt )  of 
random  phases  among  the  elemental  wavelets,  the  resultant  wavefield  clearly  obeys  Rice 
(1945)  statistics,  reducing  to  Rayleigh  statistics  for  sufficiently  strong  (in  the  aggre¬ 
gate)  scatter. 

Close  to  the  medium  and  for  sufficiently  weak  scatter,  an  identification  can  be  made 
between  the  phase  perturbations  imposed  by  the  medium  and  fluctuations  in  the  imaginary 
part  of  E.  Propagation  calculations  based  on  the  first  Born  approximation  then  yield  the 
second  moments  needed  to  characterize  the  generalized  bivariate  Gaussian  distribution  of 
x  and  y.  The  identification  of  y  with  j>,  and  of  x  with  vTT  however,  clearly  breaks  down 
for  strong  scatter,  even  close  to  the  medium. 

For  very  weak  scatter,  the  multiplicative  and  additive  views  lead  to  the  same  results 
since  fluctuations  in  E  and  in  its  logarithm  then  are  virtually  identical.  Recapping,  we 
find  (1)  that  the  multiplicative  model  is  satisfactory  for  any  degree  of  scatter  close  to 
the  medium  but  becomes  dccreasingly  satisfactory  with  increasing  distance  for  moderate  to 
strong  scatter;  (2)  that  the  additive  model  is  satisfactory  for  any  degree  of  scatter  far 
from  the  medium  but  becomes  decreasinglv  satisfactory  with  decreasing  distance  for  moderate 
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to  strong  scatter;  and  (3)  that  either  model  will  suffice  for  very  weak  scatter  at  any 
distance . 


Regimes  of  validity  for  various  signal-statistical  calculations  have  been  laid  out 
systematically  by  de  Wolf  (1969,  1975)  for  observations  inside  an  extended  medium.  His 
"saturation"  zone  corresponds  approximately  to  the  region  identified  heuristically  above 
as  not  being  well  described  by  either  the  multiplicative  or  the  additive  model  —  namely, 
the  middle-distance  regime  for  moderate  and  strong  scatter.  A  closely  related  ordering  of 
regimes  has  been  presented  by  Singleton  (1970)  for  observations  not  necessarily  within 
the  medium,  although  for  a  medium  characterized  by  a  single  scale  (Gaussian  spatial  spec¬ 
trum).  Both  these  orderings,  however,  address  only  intensity  statistics,  whereas  the 
present  work  deals  with  complex-signal  statistics.  The  latter  also  deals  directly  with 
"middle"  distances  in  which  the  Fresnel  radius  lies  between  the  inner  and  outer  scales  of  a 
power-law  scattering  medium,  both  scales  in  fact  lying  outside  the  spectral  window  filled 
by  the  data  available. 

B.  OBJECTIVE 

There  have  been  two  essential  approaches  to  obtaining  a  general  description  of  the 
signal  statistics  arising  from  scattering.  The  first  is  to  seek  an  approximate  analytical 
solution,  guided  by  experiment,  to  the  general  problem  of  defining  a  relevant  distribution 
function.  Nakagami  (1960)  has  summarized  a  remarkably  tenacious  and  quite  productive 
series  of  attacks  on  the  problem  in  Japan,  limited  however  to  description  of  intensity 
statistics.  The  second  approach  is  to  combine  the  multiplicative  acid  additive  models  in 
such  a  way  as  to  provide  a  generally  useful  hybrid  distribution.  Workers  in  the  optics 
(Wang  and  Strohbehn,  1974b) and  radio  (Fremouw,  Rino,  and  Livingston,  1976)  communities  have 
independently  proposed  very  similar  hybrid  models,  the  latter  being  somewhat  more  general 
but  requiring  more  parameters  for  its  definition. 

The  main  objective  of  the  present  work  is  to  establish  the  relative  utility  of  the 
four  leading  signal-statistical  models  —  multiplicative  (log-normal),  additive  (in  its 
generalized  Gaussian  form),  approximate  (Nakagami-m) ,  and  hybrid  (general  two-component) 

—  for  describing  the  behavior  of  radio  waves  that  have  undergone  narrow-angle  forward 
scatter  in  a  structured  plasma.  The  approach,  which  is  described  in  Section  II,  was  to 
perform  hypothesis  testing  by  means  of  chi-square  evaluation  of  distribution  functions  for 
phase  and  intensity,  employing  observations  of  the  complex  signal  received  from  a  satellite 
transmitting  through  the  ionosphere;  the  results  of  hypothesis  testing  are  presented  in 
Section  III  A.  Insight  gained  into  the  parameter  relationship  between  phase  and  intensity 
fluctuations  (scintillations)  are  described  in  Section  III  B,  and  conclusions  drawn  from 
the  work  are  presented  in  Section  IV. 


II.  EXPERIMENTAL  APPROACH 


A.  DATA  EMPLOYED 

In  May  of  1976,  satellite  P76-5  was  launched  into  a  near-circular  (1000  km),  high- 
inclination  (99.7°),  sun-synchronous  (equator  crossings  at  about  1120  and  2320,  local 
time)  orbit  with  the  DNA-002  Wideband  beacon  as  its  sole  payload.  The  Wideband  beacon 
transmits  ten  mutually  coherent  CW  signals,  ranging  from  VHF  (138  MHz)  through  UHF  (seven 
spectral  lines  centered  at  413  MHz  and  equally  spaced  11.5  MHz  apart)  to  L  band  (1239  MHz) 
and  S  band  (2891  MHz),  for  the  purpose  of  diagnosing  the  transionospheri c  radio  communi¬ 
cation  channel  and  describing  the  plasma-density  irregularities  that  perturb  it.  Obser¬ 
vations  of  dispersion  (interpreted  as  ionospheric  total  electron  content)  and  complex- 
signal  scintillation  at  the  nine  lowest  frequencies,  employing  the  S-band  signal  as  an 
essentially  unperturbed  phase  reference,  have  been  conducted  at  auroral,  mid-latitude,  and 
equatorial  stations. 

In  the  present  assessment  of  (first-order)  complex-signal  statistics,  recordings  of 
scintillating  VHF,  UHF,  and  L-band  signals  obtained  at  Poker  Flat,  Alaska  (65.4°  magnetic 
dip  latitude),  Ancon,  Peru  (0.4°  dip  latitude),  Kwajalein,  Marshall  Islands  (4.4°  dip  lati¬ 
tude),  and  Stanford,  California  (42.8°  dip  latitude)  have  been  employed.  The  quadrature 
components  of  each  signal  were  recorded  at  a  rate  of  500  samples  per  second,  after  passing 
through  low-pass  filters  having  bandwidths  of  150  Hz,  throughout  satellite  passes  of  many 
minutes  duration.  The  data  employed  in  this  analysis  were  decimated  to  100  samples  per 
second,  converted  to  intensity  and  phase,  and  detrended  to  remove  Fourier  periods  greater 
than  10  seconds  (Fremouw  at  al,  1978).  The  data  window  for  intensity  and  phase,  therefore 
corresponds  to  a  spatial  spectrum  (as  the  line  of  sight  scanned  through  the  F  layer)  be¬ 
tween  wavelengths  of  approximately  60  meters  and  30  km.  It  appears,  both  from  these  data 
and  other  experimental  results  (e.g. ,  Dyson,  McClure,  and  Hanson,  1974),  that  the  data 
window  employed  consistently  fell  between  the  inner  and  outer  scales  of  ionospheric  struc¬ 
ture. 

B.  SIGNAL-STATISTICAL  HYPOTHESES 
1.  Single-Component  Models: 

The  data  processing  employed  included  calculation  of  relevant  statistical  moments, 
computation  of  various  signal  spectra  useful  for  assessing  propagation  effects,  plots 
of  time  series  for  data-quality  control,  and  calculation  of  chi-square  fits  for  evaluation 
of  four  signal-statistical  hypotheses:  log-normal,  generalized  Gaussian,  Nakagami-m,  and 
two-component.  The  log-normal  model  is  defined  by  the  hypothesis  that  the  real  and  imag¬ 
inary  parts  of  the  complex  logarithm  of  the  received  signal  are  jointly  normal  variates. 


So-called  Gaussian  signal  statistics  follow  from  the  hypothesis  that  the  real  and 
imaginary  parts  of  the  antenna  voltage  itself,  ra ' her  than  of  its  logarithm,  are  jointly 
normal  variates.  That  is,  given 

E  =  A  e1 ''  =  x  +  iy  (5) 


then  the  joint  pdf  of  x  and  y  is 


PG  (*>y)  = 


2v,  o  ( l-o-  y 
x  y  xy 


O  z  (x-<X:)2  -2(7  a  p  (x — x-')v+0‘1  V 
y  _ x  y  xv  _  y  J 

2c  a  ‘  ( l-pz  ) 
x  v  xv 


The  mean,  variances,  and  correlation  coefficient  in  Eq.  (6)  are  analogous  to  those  in  F,q . 
(2),  but  they  are  calculated  from  the  outputs  of  a  pair  of  quadrature  detectors  rather 
than  from  the  result  of  envelope  and  phase  detection. 

in  the  general  case  of  Gaussian  statistics  represented  by  Eq .  (fS)  ,  the  individual 
pdf's  for  intensity  and  phase  are  considerably  more  complicated  than  Eqs.  (3)  and  (A). 
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Indeed,  the  intensity  pdf  is  obtainable  only  by  numerical  integration  (Fremouw  and  Rino, 
1975)  or  in  terms  of  an  infinite  series  (Beckman  and  Spizzichino,  1963),  Hatfield,  how¬ 
ever,  has  succeeded  in  obtaining  the  following  analytic  expression  for  the  phase  pdf 
(Hatfield  and  Rino,  1975): 


a  o  (1-p2  ): 

(*)  =  JL-y — sy_ 

2nD 


<x>o  a  E 
x  y 


-<X>202 


2a2o2 (1-p2  ) 
x  y  xy 


exp  <  - 


<x>2a2(l+Ea2/D) 

y  x 

2a2a2 (1-p2  ) 
x  y  xy 


-  <x>a  a  E 

x  y _ 

V2D  a  a  (1-p2  )5 
x  y  xy 


where  D  =  a2  cos2<k  -  2p  a  a  cos<f>  sin<f>  +  o2  sin2d>  (8) 

y  xy  x  y  Y  x  Y 

and  E  =  o2  cosb  -  p  o  o  sin  6,  (9) 

y  Y  xy  x  y  ' 

While  the  pdf's  for  intensity  and  phase  are  very  complicated  in  the  general  case  of 

Gaussian  statistics,  they  reduce  to  quite  familiar  forms  for  certain  special  cases.  For 

instance,  the  case  in  which  x  and  y  have  equal  variances  (a2  and  )  and  are  uncorrelated 

(pxy  =  0)  produces  the  well-known  Rice  statistics  obeyed  by  a  steady  signal  accompanied  by 

noise.  When  the  noise  totally  dominates  the  signal  (a2  +  a2  =  <T>),  we  have  Rayleigh  sta- 

x  y 

tistics. 

A  more  general  approximation  to  Gaussian  statistics  has  been  worked  out  by  Nakagami 
(1960)  by  means  of  a  semi-empirical  procedure.  The  so-called  Nakagami-m  distribution  for 
intensity,  which  approaches  the  Rice  distribution  as  m  approaches  unity  from  above  and 
equals  the  Rayleigh  distribution  when  m  =  1,  is 


m  an—  1 
m  r 


Pv  (I)  =  - - ~ 

r(m)<I> 


-ml /<!> 


;I2>-<I>2 


The  Nakagami  distribution  has  attained  considerable  favor  among  observers  of  intensity 
scintillation  (Whitney,  1974),  since  it  is  independent  of  any  phase-related  parameter. 

By  the  same  token,  there  is  no  phase  pdf  accompanying  the  Nakagami  distribution,  nor  is 
there  a  related  description  of  the  joint  statistics  of  intensity  and  phase. 

2.  Two-Component  Model: 

Even  the  generalized  Gaussian  statistics  yield  only  a  limited  description  of  phase. 


signal-statistical  models.  The  two-component  model  tested  in  the  present  work  is  that  put 
forth  by  Fremouw,  Rino,  and  Livingston  (1976),  in  which  the  total  scintillating  signal,  E, 
is  modeled  as  the  product  of  a  dif fractively  scattered  component,  Eg  ,  and  a  refractively 
focused  component,  E^ ,  such  that 

E  =  A  e1^  =  EgEf  =  (xs+iys  )exp(xf+iij>f  ) .  (12) 

In  this  two-component  model,  it  is  postulated  that  Eg  obeys  generalized  Gaussian  sta¬ 
tistics  [Eq.  (6)]  and  that  E^  obeys  log-normal  statistics  [Eq.  (2)].  A  heuristic  rationale 
for  the  two-component  model  can  be  made  by  considering  a  wave  that  passes  through  a  region 
of  large-scale,  phase-perturbing  structure  and  then  through  a  layer  of  smaller-scale  scat- 
terers.  The  phase  of  the  wave  would  have  undergone  a  random  walk  through  the  first  region, 
and  weak  focuses  and  defocuses  would  develop  in  subsequent  propagation.  It  is  postulated 
that  this  process  leads  to  log-normal  statistics. 

Upon  passing  through  the  second  layer,  the  wave  would  emerge  with  a  scattered  (co¬ 
herent)  and  a  non-scattered  (noncoherent)  part.  The  latter  would  retain  log-normal  statis¬ 
tics,  but  the  former  would  be  disrupted  by  the  diffraction  process,  leading  to  Gaussian 
statistics.  A  tacit  assumption  in  the  two-component  model  is  that  the  resulting  statistics 
are  independent  of  whether  the  small-scale  scatter  takes  place  before,  after,  or  concur¬ 
rently  with  the  large-scale  phase  perturbations.  The  two  components  are  taken  to  be  mul¬ 
tiplicative  because  signal  components  behaving  qualitatively  as  those  postulated  in  the 
model  were  isolated  by  means  of  a  coherent-AGC  process  which  is  tantamount  to  dividing  E 
by  Ej ,  as  illustrated  by  Fremouw  et  al  (1978). 

A  basic  postulate  in  the  two-component  model  is  that  Eg  and  E^  are  statistically  in¬ 
dependent.  Experimentally,  this  is  ensured  by  isolating  the  two  components  by  means  of  a 
very  sharp  filter  (a  ten-pole  Butterworth  in  the  present  work).  In  terms  of  propagation 
theory,  this  experimental  step  corresponds  to  a  partitioning  of  the  scattering  medium  in 
such  a  way  that  the  scale  for  changes  in  the  transverse  gradient  of  in  E^.  are  much  larger 
than  the  scale  for  changes  in  the  transverse  gradient  of  Eg .  Under  such  a  partitioning, 
the  parabolic  wave  equation  can  be  solved  independently  for  the  two  components,  as  has 
been  carried  out  by  Rino  (Fremouw  and  Rino,  1978). 

Under  the  postulate  of  statistical  independence,  it  is  straightforward  to  show  that 


<EE*>  =  (<x  >2  +  a2)  exp  (<xf:>  +  a2  )  =  1 
s  s  T  Xj 

and  <EE>  =  (<xg  >2  +  Bg  )  exp  [2<x^>  +  ] 

A  o  o 

where  B_  =  az  -of  +  2ip  a  a 


f  X 


Af  vf  Xf't’f  Xf  4>f 

a|  =  a2+a2  (17) 

s  ^  s 

and  B  =  a2  -  a2  +  2ip  a  a  (18) 

sxy  xyxy 

s  s  s  s  s  s 

and  where  we  have  used  the  fact  that  the  experimental  separation  of  Es  and  E^  results 
approximately  in  <ys  >  =  <4>^  >  =  0  and  where  conservation  of  energy  is  invoked  as  the  nor¬ 
malization  condition  [Eq.  (14)].  It  also  follows  that 


<&*->  —  <<b~>  + 


and  S2  =  t  (S?  )  +1  ]  [(S2),+l]  -  1 


where  (S^ )  ^  are  the  component  intensity  scintillation  indices  [defined  according  to  Eq. 

(11)]. 

Cursory  inspection  of  Eq.  (20)  might  suggest  that  the  two-component  model  predicts 
occurrence  of  values  considerably  in  excess  of  unity  (for  instance,  =\  3  for  = 

S4f  =  1).  While  the  model  would  admit  such  a  situation,  it  does  not  in  fact  predict  it 
since  it  does  not  address  the  question  of  a  parameter  dependence  between  and  .  The 
largest  value  of  yet  observed  in  the  Wideband  experiment  is  1.2,  and  detailed  inspection 
of  signal  behavior  suggests  that  only  the  coherent  part  of  Eg  may  be  susceptible  to  signif¬ 
icant  focusing  and  defocusing,  which  would  put  a  limit  on  simultaneously  large  (near  unity) 

values  of  S.  and  S.  ,  . 

4s  4f 

For  purposes  of  hypothesis  testing,  we  need  the  pdf's  for  phase  and  intensity  under 
the  two-component  model.  For  phase,  the  problem  is  trivial  because  =  <)>£+^s,  where 
'tig  =  tan  ^(ys/xs)  and  where  and  i(>  are  statistically  independent.  Thus,  the  pdf  for  the 
composite  phase  is  the  following  convolution: 

Pc$U)  =  -A  (*s)pn(*-*s)d*3  (21) 

where  is  the  Hatfield  distribution  [ Eq .  (7)]  and  p^  is  the  normal  distribution  [Eq.  (4H 


I 
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Similarly,  from  the  statistical  independence  of  ^  and  Xg (  =  i  in  I  ) ,  one  can  show, 
by  means  of  a  series  of  coordinate  transformations,  that  the  pdf  for  the  composite  inten¬ 
sity  in  the  two-component  model  is 

Pd  (I)  =  -/C\l(Is)pfc(f)dIs  (22) 

where  I  =  x2  +  y2  and  where  p_T  is  the  intensity  pdf  arising  from  numerical  integration  of 

S  S  S  (jJ. 

Eq.  (6)  and  p  is  the  log-normal  pdf  defined  in  Eq.  (3). 

More  than  just  the  individual  pdf's,  however,  the  two-component  hypothesis  provides  a 
joint  distribution  for  intensity  and  phase.  Employing  reasoning  similar  to  that  leading  to 
Eqs.  (21)  and  (22)  and  invoking  statistical  independence  of  the  two  components,  one  finds 
the  following  joint  pdf: 

2  X  v 

Pc  =  27_JTe  Pg  (e  Scos  <f>g,e  sin  ^s)Pn  (i^nl-)^  ,<(i-<(.s  jd^d^  (23) 

which  is  essentially  a  two-dimensional  convolution  of  the  two  bivariate  normal  distrib¬ 
utions  defined  in  Eqs.  (2)  and  (6). 

C.  DATA  PROCESSING 
1.  Sample  Selection: 

To  obtain  a  representative  set  of  data  samples  for  detailed  signal-statistical  anal¬ 
ysis,  the  entire  collection  of  (singly)  detrended  VHF  strip  charts  (cf.  Figure  6  of  Fremouw 
et  al,  1978)  obtained  during  approximately  the  first  year  of  the  Wideband  experiment  was 
inspected.  The  inspection  was  aimed  at  obtaining  about  one  hundred  data  samples  of  at 
least  20  seconds  duration  during  which  the  received  signals  displayed  reasonable  sta¬ 
tistical  stationarity .  Selection  criteria  included  an  attempt  to  find  samples  reasonably 
well  distributed  in  satellite  elevation  ingle  and  scintillation  disturbance  level  and  to 
have  representation  from  all  four  receiving  stations:  Ancon,  Kwajalein,  Poker  Flat,  and 
Stanford. 

Table  1  shows  the  distribution  of  record  segments  selected  (ranging  in  duration  from 
20  to  85  seconds),  for  all  of  wiiich  VHF  data  were  processed.  For  the  more  disturbed  seg¬ 
ments,  meaningful  analysis  could  be  conducted  at  higher  frequencies,  and  the  number  of  VHF, 
UHF,  and  L-band  data  segments  analyzed  from  each  station  is  given  in  Table  2. 
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Table  1.  Distribution  of  VHF  dual-detrend  data  sets  processed 


Location 

Elevation 

Weak 

Disturbance  level 

Moderate 

Strong 

Low 

(  <20°) 

1 

2 

Ancon 

Medium 

(20-50°) 

1 

1 

4 

High 

(  >50°) 

3 

Low 

(  <20°) 

2 

3 

5 

Kwajalein 

Medium 

(20-50°) 

4 

5 

3 

High 

(  >50°) 

1 

Low 

(  <20°) 

1 

1 

1 

Poker  Flat 

Medium 

(20-50°) 

2 

5 

2 

High 

(  >50°) 

2 

1 

Low 

(  <20°) 

1 

3 

3 

Stanford 

Medium 

(20-50°) 

1 

3 

High 

(  >50°) 

1 

Table  2.  Number 

of  dual-detrend 

data  sets  processed 

Location 

VHF 

UHF 

L-band 

Total 

Ancon 

12 

9 

6 

27 

Kwajalein 

23 

11 

4 

38 

Poker  Flat 

15 

6 

4 

25 

Stanford 

12 

4 

16 

62 

30 

14 

106 

14 


2. 


The  Dual-Detrend  Processor: 

Under  this  contract,  the  dual  detrending  procedure  developed  by  Fremouw  and  Living¬ 
ston  (1976)  was  employed  for  analysis  of  che  data  sets  summarized  in  Tables  1  and  2,  in 
order  to  test  the  two-component  hypothesis  and  the  three  single-component  hypotheses.  The 
processor  also  provides  a  wealth  of  complementary  information  on  signal  behavior,  per¬ 
forming  numerous  computational  functions  and  delivering  several  numerical  and  graphical 
outputs . 

Central  to  the  processing  performed  is  computation  of  a  number  of  first-order,  signal- 

statistical  moments.  For  the  composite  (ie,  0.1-Hz  singly  detrended)  signal,  the  following 

moments  are  computed:  the  mean  intensity, <T>;  the  mean  log-amplitude,  <x>;  the  mean  real 

and  imaginary  parts,  <x>  and  <y>  (=0);  the  normalized  standard  deviation  of  intensity,  ; 

the  standard  deviation  of  log-amplitude,  a  ;  the  standard  deviation  of  phase,  o, ;  the 

X  Y 

variances  of  the  real  and  imaginary  parts,  o2  and  a2;  the  total  signal  variance,  a2 

x  y  t 

=  a2  +  a2;  and  various  combinations  of  these  moments, 
x  y 

For  use  in  testing  the  two-component  model,  several  analagous  quantities  are  computed 
for  each  of  the  two  components,  which  are  separated  by  independently  passing  the  composite 
phase  and  log-amplitude  through  a  low-pass,  ten-pole  Butterworth  filter  having  a  3-dB 
cutoff  at  0.4  Hz,  storing  the  results,  and  then  subtracting  them  from  their  nonfiltered 
counterparts.  The  smoothed  results  comprise 


Ef  =  exp  (Xj+iij’f)  =  exp  (£+i$) 


(24) 


where  “  denotes  the  smoothing  process,  and  the  residuals  comprise 


Eg  =  xs+iyg  -  exp  [(x-8l+ i  (<>-<>)]  . 


(25) 


The  component  moments  computed  include  the  following  "six  sigmas"  needed  to  define  the 


component  pdf's:  o,<f,  n^f,  «Xf  +  f  =VP> 


J  \[  f AP  ^  ^  ^  j  U  \ .  p  O  U  • 

Ar<J>r  Xt  <J>t  xs  y s  xsys  v  xsys  xs  ys 
addition,  quantities  derived  from  the  six  sigmas,  such  as  those  defined  in  Eq .  (16),  (17), 

and  (18),  are  calculated  for  comparison  of  signal  behavior  with  that  predicted  by  the  two- 
component  model,  such  as  that  described  by  Eqs.  (13),  (14),  and  (15). 

Second-order  moments  of  the  composite  signal,  in  the  form  of  several  spectra  and 
autocorrelation  functions,  are  computed  and  output  graphically.  Specifically,  plots  are 
obtained  of  the  power  spectra  and  autocorrelation  functions  of  intensity,  I  and  phase,  <f, 
and  of  the  two-siued  (positive  and  negative  frequencies)  Doppler-spread  power  spectrum  of 
the  received  signal,  E.  Since  the  signal  fluctuations  recorded  arise  from  scanning  of  a 
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complex  diffraction  pattern  across  the  receiving  antenna,  the  Doppler  spectrum  also  may 
be  viewed  as  a  one-dimensional  cut  through  the  two-dimensional  angular  spectrum  of  plane 
waves  arriving  at  the  ground  (Booker,  Ratcliffe,  and  Shinn,  1950). 

Hypothesis  testing  is  carried  out  by  forming  histograms  of  intensity  and  phase,  sep¬ 
arately,  from  the  (composite)  complex  signals  recorded  and  then  comparing  the  histograms 
with  pdf's  computed  on  the  basis  of  the  four  signal-statistical  hypotheses.  The  theoretical 
pdf's  are  calculated  from  the  relevant  moments  computed  from  each  data  set.  For  testing  the 
log-normal  hypothesis,  <x>(*0)  and  are  inserted  in  Eq.  (3),  and  is  employed  in  Eq. 

(4).  For  testing  the  generalized  Gaussian  hypothesis,  <x>,  o^,  o  ,  and  are  used  in 

Eqs.  (7),  (8),  and  (9),  and  Eq.  (6)  is  integrated  numerically  to  obtain  p  (I)  from  the 

Gi 

same  moments.  (Note  that  the  normalization  condition  reduces  the  number  of  parameters  to 

three,  since  <x>2  =  o2+cr2.) 

x  y 

The  Nakagami  distribution  is  calculated  from  Eqs.  (10)  and  (11),  employing  the  ob¬ 
served  intensity's  mean,  <I>  (=1  by  normalization)  and  variance,  <I2>,  For  testing  the 
two-component  model,  the  six  sigmas  defined  after  Eq.  (25)  are  employed  to  calculate  the 
component  distributions  appearing  in  Eqs.  (21)  and  (22),  and  the  convolutions  are  carried 
out  numerically. 

For  visual  comparison,  the  intensity  and  phase  histograms  and  the  theoretical  pdf's 
are  output  graphically.  Quantitative  hypothesis  testing,  however,  is  performed  by  calcu¬ 
lating  the  chi-square  difference  between  each  histogram  and  each  corresponding  theoretical 
pdf.  An  additional  graphical  output  gives  a  comparison  of  the  histograms  of  X^,  4^,  xg, 
and  y^  with  normal  curves  computed  from  the  relevant  means  and  standard  deviations.  No 
quantitative  use  is  made  however,  of  this  information  about  the  components  of  the  two- 
component  model. 
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III.  RESULTS 


A.  HYPOTHESIS  TESTS 

1.  Chi-Square  Values 

Complete  dual-detrend  processing,  including  chi-square  testing  of  the  four  signal- 
statistical  hypotheses,  was  conducted  on  101  of  the  106  data  samples  enumerated  in  Table  2. 
(The  other  five  were  processed  as  test  cases  during  software  development,  before  the  chi- 
square-test  routine  was  operational.)  With  a  few  individual  exceptions  (mainly  involving 
software  limitations  in  generating  the  Nakagami-m  pdf  for  extremely  small  S^  ) ,  valid  chi- 
square  tests  were  obtained  for  all  distributions.  In  this  section,  we  present  and  discuss 
the  behavior  of  x2  as  an  indicator  of  the  efficacy  of  the  four  signal-statistical  models. 

While  it  is  to  be  hoped  that  a  single  model  can  be  employed  for  all  observing  condi¬ 
tions,  one  might  expect  to  find  differences.  Accordingly,  the  test  results  have  been 
grouped  for  comparison  by  observing  station  and  frequency.  One  might  also  expect  geomet¬ 
rical  differences  (e.g.,  a  trend  with  elevation  angle),  but  the  geometrical  situation  is 
sufficiently  complicated  that  a  simple  sorting  of  x2  probably  would  not  be  productive. 
Geometrical  effects  do  play  an  important  part  in  establishing  the  moments  (e.g.,  the  phase 
and  intensity  scintillation  indices,  and  )  by  which  the  signal  statistics  are  para¬ 
meterized,  and  this  effect  will  be  explored  in  Section  III  B. 

Table  3  gives  the  mean  values  and  standard  deviations  of  x2  for  the  various  models  for 
each  of  the  four  stations  and  three  observing  frequencies.  (VHF  is  138  MHz;  UHF  includes 
measurements  at  both  379  and  413  MHz;  and  L-band  is  1239  Mhz.)  The  numbers  in  square 
brackets  preceding  each  mean  value  indicate  the  ranking  of  the  model  in  that  column  for  the 
data  sets  in  that  row.  (I.e.,  [1]  denotes  the  smallest  mean  x2»  [2J  the  next  smallest,  and 
so  forth.) 

Close  examination  of  Table  3  reveals  several  aspects  of  the  signal-statistical  behavior 
of  scintillating  complex  signals.  Let  us  peruse  it  first  for  trends  with  frequency.  Taking 
note  of  the  model  rankings  for  phase,  we  find  that  for  no  station  does  the  ranking  change 
in  going  from  VHF  to  UHF,  nor  does  it  change  at  Ancon  in  going  to  L-band.  The  differences 
observed  between  the  various  mean  x‘  values  for  L-band  at  Kwajalein  are  very  small,  com¬ 
pared  with  the  corresponding  standard  deviations,  and  those  at  Poker  Flat  are  only  slightly 
less  so.  (Furthermore,  the  L-band  data  segments  from  Kwajalein  and  Poker  Flat  represent 
only  17.  of  the  total  population.) 

We  conclude  that  there  is  no  statistically  significant  change  in  the  ranking  of  phase- 
statistical  models  with  changing  frequency  (and,  therefore,  with  changing  Fresnel-zone 
size).  As  had  been  inferred  earlier  (Fremouw  et  al,  1978)  from  the  frequency  dependence  of 

o  ,  which  is  essentially  independent  of  Fresnel-zone  size  in  most  rases,  the  first-order 
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Table  3.  Means  and  standard  deviations  of  x2 
([n]  denotes  ranking) 


Frequency 

(No.  of  Log-Normal 

X2  Values)  _ 


Generalized 

Gaussian 


Two-Component  Nakagami-m 


VHF  (12) 


Phase 

til 

5.41+6.46 

[33 

25.06+12.53 

[2] 

6.2916.77 

Intensity 

[4] 

104.46+99.20 

[2] 

7.87+5.83 

[3] 

14.38110.96 

[1] 

7.11+5.36 

UHF  (9) 

Phase 

[1] 

3.26+1.69 

[3] 

10. 02±9.82 

[2] 

3.68+2.07 

Intensity 

[3] 

39.11±31.88 

[4] 

45. 19±53. 16 

[11 

8.90+5.59 

[21 

9.13+6.18 

L-band  (6) 

Phase 

[1] 

3.43+1.82 

[3] 

4.87+2.24 

[2] 

3.93+1.99 

Intensity 

[1] 

5.30+2.55 

[4] 

96.62H15.42 

[31 

10.55+13.24 

[2] 

8.06+4.52 

VHF  (23) 

Phase 

[1] 

6.24+6.72 

[3] 

10.1519.34 

[2] 

7.15+7.65 

Intensity 

[4] 

24.45±27.08 

[11 

6.67+7.38 

[3] 

10. 92±8.89 

[2] 

7.27+7.60 

UHF  (11) 

Phase 

[1] 

4.16+4.48 

[3] 

5.59+6.65 

[2] 

4.60+5.08 

Intensity 

[4] 

15.58+21.07 

[2] 

6.7518.31 

[33 

10.35+6.44 

U1 

6.46+7.25 

L-band  (3) 

Phase 

[2] 

1.86+1.32 

[3] 

1.95+1.25 

[11 

1.7811.27 

Intensity 

[2] 

3.20±1.69 

[4] 

6.76+7.11 

[3] 

3.28+1.42 

[11 

2.76 

(sincle  value) 

VHF  (12) 

Phase 

[1] 

2.38±1.27 

[3] 

8.47+5.50 

[2] 

2.73+1.31 

Intensity 

[2] 

20.86+22.00 

[4] 

149.02+159.32 

[3] 

29. 88156.13 

[13 

8. 95_  8.92 

UHF  (5) 

Phase 

[11 

2.38+1.42 

[3] 

7.1614.19 

[2] 

2.4311.43 

Intensity 

[2] 

20.97+18.91 

[4] 

288.611162.00 

[3] 

59.11+35.49 

HI 

20.21119.56 

L-band  (4) 

Phase 

[2] 

1.7210.74 

[13 

1.51+0.42 

[33 

1.80+0.61 

Intensity 

[2] 

10.03+5.79 

[4] 

262.97+181.09 

[31 

14.75+9.14 

[  1?  ](sof  twari' 

limit) 

VHF  (12) 

Phase 

U3 

4.49+5.21 

[3] 

9.1815.70 

[21 

4.91+5.72 

Intensity 

[3] 

14.59+9.46 

[4] 

174.42+141.66 

[2] 

12.04+13.51 

[1] 

8.88+6.43 

VHF  (4) 

Phase 

[11 

3.79+3.63 

[3] 

4.14+4.06 

[2] 

3.84+3.16 

Intensity 

[31 

8.49+6.01 

[4] 

40.35+36.43 

[2] 

6.81+4.24 

[11 

6.66+5.13 

statistics  of  phase  are  quite  durable  in  post-scattering  propagation.  There  is  evidence, 
however,  of  measurable  change  in  the  second-order  statistics  —  namely,  some  decrease  in 
the  power-law  index  of  the  phase  spectrum  as  intensity  scintillation  develops. 

Regarding  intensity  statistics,  we  find  no  differences  in  ranking  when  changing  fre¬ 
quency  at  Poker  Flat  or  Stanford.  At  the  two  equatorial  stations,  where  there  were  more 
cases  of  moderate  to  strong  scintillation  extending  to  higher  frequencies,  however,  one 
can  discern  a  consistent  frequency  trend  for  test  results  from  the  two  classical  models. 
That  is,  there  is  a  trend  for  the  log-normal  distribution  to  improve  its  ranking  and  for 
generalized  Gaussian  statistics  to  deteriorate  with  increasing  frequency  (decreasing 
Fresnel-zone  size).  This  behavior  is  what  might  be  expected  for  transmission  through  a 
remote  scattering  medium  having  an  outer-scale  cutoff  in  its  spatial  spectrum,  and  it  may 
be  understood  in  terms  of  imposition  of  an  effective  outer  scale  during  separation  of 
scintillations  from  dispersive-phase  trends. 

A  notable  feature  in  Table  3  is  that  the  relative  efficacy  of  the  two  classical 
models  is  different  at  different  stations.  In  particular,  better  intensity  fits  are  ob¬ 
tained  at  VHF  with  the  generalized  Caussian  hypothesis  than  with  the  log-normal  hypothesis 
for  the  two  equatorial  stations  (Ancon  and  Kwajalein),  but  the  reverse  is  true  for  the  two 
higher-latitude  stations.  This,  too,  may  be  understood  in  terms  of  an  imposed  outer 
scale,  which  turns  out  to  be  geometry-dependent,  as  will  become  clear  in  Section  III  B. 

A  key  point  to  be  noted  is  that  the  two-component  model  represents  an  improvement 
over  the  less  favored  of  the  two  classical  models  for  every  combination  of  frequency  and 
observing  station  (eleven  independent  data  populations)  for  both  intensity  and  phase. 
Moreover,  the  two-component  model  provides  mean  x‘:  values  that  are  nearly  as  favorable  as 
(in  a  few  instances,  better  than)  the  values  for  the  better  of  the  two  classical  hypo¬ 
theses,  when  measured  in  terms  of  the  relevant  standard  deviation  of  x2  or  tine  difference 
in  mean  x2  between  the  two  classical  models.  Thus,  the  two-component  hypothesis  provides 
a  generally  useful  model  for  the  complex-signal  statistics  of  scintillation,  including  the 
joint  pdf  of  intensity  and  phase  [Eq.  (23)]. 

The  other  key  point  is  the  consistency  with  which  the  Nakagami-m  distribution  pro¬ 
vides  a  good  fit  to  intensity  data.  In  terms  of  mean  x2»  the  Nakagami  distribution  pro¬ 
vided  the  best  description  of  intensity  statistics  in  seven  of  the  ten  Table-3  categories 
for  which  it  was  reliably  tested,  and  it  was  second  best  of  the  four  hypotheses  in  the 
remaining  three  categories.  The  efficacy  of  the  Nakagami-m  distribution  is  illustrated 
also  in  Table  4,  which  summarizes  model  rankings  for  both  phase  and  intensity,  giving  the 
total  number  of  times  that  eacli  model  achieved  each  rank.  Nakagami-m  provided  the  smal¬ 
lest  value  for  intensity  in  32  cases  out  of  the  K3  in  which  all  four  models  were  re- 
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liably  tested,  ranked  second  in  38  cases,  was  third  in  ten  cases,  and  provided  the 
least  favorable  fit  in  only  three  tests. 


Table  4.  Model  Rankings 
(Number  of  cases) 


Rank 

Log-normal 

Generalized 

Gaussian 

Two-Component 

Nakagami-m 

1 

Phase 

54 

23 

22 

- 

Intensity 

11 

28 

12 

32 

2 

Phase 

Intensity 

38 

10 

16 

11 

45 

24 

38 

i 

3 

Phase 

Intensity 

7 

22 

60 

8 

32 

43 

10 

4 

Phase 

Intensity 

40 

36 

4 

3 

comparison,  the  same 

83  tests  gave 

first  rank  to 

the  two-component 

model  12  times. 

the  log-normal  model  11  times,  and  to  the  generalized  Gaussian  model  28  times;  the  latter 
ranked  second  11  times,  third  8  times,  and  last  36  times.  We  conclude  that  for  intensity 
only,  Nakagami  provides  a  more  consistently  good  statistical  distribution  than  either  log¬ 
normal,  generalized  Gaussian  (when  the  theoretical  distribution  is  calculated  from  measured 
moments  of  the  real  and  imaginary  signal  components),  or  the  two-component  hybrid  model  of 
Fremouw,  Rino,  and  Livingston  (1976). 

Rino,  Livingston,  and  Whitney  (1976)  concluded  from  analysis  of  intensity-only  data 
that  the  statistics  underlying  complex-signal  scintillation  are  generalized  Gaussian  but 
highly  non-Rician.  They  obtained  better  fits  to  the  generalized  Gaussian  model  than  to  the 
log-normal  model,  but  their  method  of  testing  the  Gaussian  hypothesis  allowed  adjustment  of 
three  free  parameters,  whereas  the  log-normal  distribution  was  totally  specified  by  a  mea¬ 
sured  moment;  therefore,  their  conclusion  in  favor  of  generalized  Gaussian  statistics  was 
not  a  strong  one.  They  did,  nonetheless,  use  their  best— fit  parameters  to  infer  that  scin¬ 
tillation  is  usually  accompanied  by  greater  variance  in  the  imaginary  signal  component  than 
in  the  real  component.  The  Wideband  experiment  has  provided  direct  evidence  in  support  of 
the  latter  conjecture,  but  it  has  demonstrated  inadequacy  of  the  generalized  Gaussian  model 
to  provide  a  full  description  of  complex-signal  statistics. 

The  Nakagami  distribution  was  largely  ignored  by  Rino,  Livingston,  and  Whitney  (1976) 
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partly  because  they  viewed  it  as  virtually  identical  to  the  Rice  distribution.  Numerical 
calculations  by  the  present  authors,  however,  have  disclosed  differences  between  the  Rice 
and  Nakagami  distributions  for  intermediate  values  of  (0.4  to  0.8,  say)  at  least  as 
large  as  the  differences  between  the  log-normal  and  generalized  Gaussian  distributions  com¬ 
pared  by  the  former  authors.  Indeed,  they  might  have  found  Nakagami  as  satisfactory  as 
generalized  Gaussian  without  having  to  vary  free  parameters  to  obtain  good  fits. 

The  shortcoming  of  the  Nakagami  model,  of  course,  is  that  the  m  distribution  has  no 
counterpart  for  phase.  Referring  again  to  Table  3,  we  see  that  the  log-normal  hypothesis 
provided  the  smallest  mean  x2  value  for  phase  in  all  11  data  categories  contained  in  the 
table.  Table  4  also  shows  that  phase  usually  is  normally  distributed.  Thus,  a  generally 
useful  model  for  complex-signal  statistics  would  be  one  that  reduces  to  the  product  of  the 
normal  distribution  for  phase  and  the  Nakagami  distribution  for  intensity  when  correlation 
between  $  and  I  is  zero  but  that  accounts,  in  general,  for  the  correlation  that  is  observed 
experimentally  (Fremouw,  Rino,  and  Livingston,  1976;  Fremouw  et  al,  1978).  Until  such  a 
bivariate  distribution  is  developed  and  tested,  the  two-component  model  can  serve  to  de¬ 
scribe  the  full  statistics  of  scintillating  complex  signals. 

A  final  point  about  intensity  distributions  is  illustrated  in  Figure  1,  which  is  a 
plot  of  x2  values  against  the  intensity-scintillation  index,  ,  observed  at  Kwajalein  (all 
frequencies).  In  general,  the  raw  x2  values  increase  somewhat  with  increasing  index,  which 
is  a  measure  of  the  spread  in  the  histogram  of  observed  intensity  values.  The  significant 
point  is  that  the  increase  in  x2  with  increasing  S,  is  much  more  pronounced  for  the  log- 
normal  distribution  than  for  the  other  three.  This  is  a  clear  indication  that  the  log¬ 
normal  distribution,  which  arises  from  the  multiplicative  propagation  model,  is  less  effi¬ 
cacious  in  the  "saturation"  regime  than  are  distributions  arising  from  the  additive  model, 
which  is  the  basis  of  the  Nakagami  and  generalized  Gaussian  hypotheses  and  which  plays  a 
dominant  role  in  development  of  the  intensity  statis  ics  described  by  the  two-component 
hypothesis.  This  effect  may  be  understood  in  terms  of  a  randomizing  of  phases  in  the 
broader  angular  spectrum  that  accompanies  stronger  scatter. 

2.  A  Representative  Example 

Section  III  A  1  deals  with  the  overall  statistical  behavior  of  approximately  100  record 
segments.  Aside  from  the  trends  identified  in  that  section,  there  is  a  good  deal  of  var¬ 
iation  from  data  set  to  data  set,  as  is  indicated  by  the  large  standard  deviations  of  x2 
enumerated  in  Table  3.  In  this  section,  we  describe  the  behavior  of  a  single  80-second 
data  segment  as  a  means  of  illustrating  salient  features  in  the  larger  data  population. 

The  segment  chosen  is  from  a  pass  over  Ancon,  starting  at  0420:19  GMT  on  4  March  1977 
(2320:19  local  civil  time  on  3  March)  when  the  satellite  was  at  an  elevation  of  about  70° 
at  an  azimuth  near  170°.  The  intensity  scintillation  index,  S,  ,  was  0.90  (m  =  1.05)  at  VHF, 
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0.66  (m  -  1.23)  at  UHF ,  and  0.13  (m  -  2.77)  at  L  band. 

Figure  2,  which  presents  time  plots  of  the  phase,  intensity,  and  quadrature  compo¬ 
nents  of  the  VHF  signal,  E,  plus  a  scatter  plot  on  the  complex  plane,  suggests  that  the 
ionosphere  was  behaving  essentially  as  a  Rayleigh  communication  channel  at  VHF.  Computed 
moments  reveal  that  a  little  over  1%  of  the  received  power  was  in  a  steady  (over  ten  sec¬ 
onds)  component  and  that  slightly  over  half  (52%)  of  the  fluctuation  power  resided  in  the 
imaginary  component. 

Figure  3  shows  the  time-series  and  complex-plane  behavior  of  the  two-component  model's 
scatter  component,  Eg  ,  which  contained  a  coherent  part  bearing  88%  of  the  power  and  a  non¬ 
coherent  part  whose  power  was  56%  imaginary  and  44%  real.  The  two-component  model's  focus 
component,  E^  ,  is  illustrated  in  Figure  4.  The  standard  deviation,  ,  of  the  imaginary 
part  of  its  complex  logarithm  was  2.83  radian,  which  accounted  for  most  of  the  composite 
phase  fluctuation,  »  3.32  radian,  and  it  imposed  an  intensity  modulation  such  that 
a  ,  =  0.20. 

xf 

The  results  of  hypothesis  testing  for  this  case  are  shown  in  Figure  5,  which  displays 
the  observed  intensity  and  phase  histograms  and  the  theoretical  pdf's  calculated  from  the 
measured  moments.  In  terms  of  ranking  on  the  basis  of  raw  x2  values,  the  behavior  Is 
typical  of  VHF  behavior  at  Ancon.  (Compare  Table  3.)  The  figure  illustrates  well  the 
failings  of  the  log-normal  distribution  for  intensity  in  the  saturation  regime:  namely, 
that  it  is  too  peaked  and  underestimates  the  probability  density  for  deep  fades. 

The  generalized  Gaussian  hypothesis  produces,  for  this  event,  essentially  a  Rayleigh 
distribution,  which  does  not  underestimate  deep  fades.  Indeed,  when  is  computed  from 
the  measured  moments  of  x  and  y,  it  consistently  overestimates  deep  fades.  This  defic¬ 
iency  worsens  as  the  Fresnel-zone  size  decreases,  as  may  be  seen  in  Figure  6,  which  shows 
histogram  fits  for  the  UHF  signal.  The  deficiency  may  be  understood  from  Figure  7,  which 

displays  complex-plane  scatter  plots  for  E,  E  ,  and  E,  at  VHF,  UHF,  and  L  band.  As  the 

s  t 

Fresnel  radius  decreases  (increasing  frequency),  the  measured  value  of  stems  decreas- 
ingly  from  intensity  scintillations  and  increasingly  from  phase  fluctuation  (in  relative 
terms).  The  generalized  Gaussian  model  cannot  accommodate  such  behavior. 

Figure  8,  which  shows  the  fits  for  L  band,  demonstrates  the  tendency  for  all  the  pdf's 
to  become  identical  for  weak  scatter.  In  the  limit  of  very  weak  scatter,  all  the  models 
become  satisfactory;  at  L  band  in  the  present  case,  however,  there  is  still  sufficient 
phase  fluctuation  that  the  measured  value  of  produces  a  poor  generalized-Gaussian  inten¬ 
sity  fit. 

The  fact  that  good  fits  can  be  made  to  highly  selected  intensity  data  by  optimizing 
free  parameters  in  the  generalized-Gaussian  model,  as  was  done  by  Rino,  Livingston,  and 
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PROBABILITY  DENSITY  -  relative 


Histograms  of  phase  and  intensity  of  the  VHF  signal  shown  in  Figure  2 ,  com¬ 
pared  with  probability  density  functions  (pdf's)  calculated  on  the  basis  of 
four  signal-statistical  hypotheses. 


Figure  6 


Histograms  of  phase  and  intensity  of  the  UHF  (379-MHz)  signal  received  simul¬ 
taneously  with  the  VHF  signal  shown  in  Figure  2,  compared  with  calculated 
pdf ' s. 
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Figure 


7. 


Complex-plane  scatter  plots  for  the  VHF  (top),  UHF  (middle),  and  L-Band  (bottom) 
signals  received  at  Kwajalein  during  the  80  seconds  beginning  at  0420:19  GMT  on 
4  March  1977.  T.eft:  composite  signil,  E.  Center:  focus  component,  E^. 

Right:  scatter  component,  Eg. 
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PROBABILITY  DENSITY  -  relative 


Whitney  (1976)  suggests  that  a  component  of  complex-signal  scintillation  might  obey  such 
statistics.  This  conjecture  formed  part  of  the  experimental  basis  for  the  two-component 
hypothesis  of  Fremouw,  Rino,  and  Livingston  (1976). 

Figure  5  illustrates  the  fact  that,  when  the  intensity  distribution  approximates  pGI , 
the  two-component  model  retains  the  essential  characteristics  of  generalized-Gaussian  sta¬ 
tistics  as  they  pertain  to  intensity,  providing  a  fit  that  is  better  than  p^  and  nearly  as 
good  as  pGj .  Figures  6  and  8  illustrate  the  durability  of  the  two-component  model  as  the 
Fresnel-zone  radius  decreases  toward  a  regime  where  log-normal  statistics  are  satisfactory. 

Figure  9  compares  histograms  of  xg  and  yg  and  of  and  ^ ,  for  L  band,  with  normal 
distributions.  Its  top  two  plots,  together  with  the  complex-plane  scatter  diagram  for  Eg  at 
L  band  in  Figure  7,  reveal  a  truly  generalized-Gaussian  component.  It  is  because  this  com¬ 
ponent  usually  dominates  intensity  scintillations  that  Rino,  Livingston,  and  Whitney  (1976) 
were  able  to  obtain  good  intensity  fits  on  the  basis  of  the  generalized-Gaussian  hypothesis. 
Isolation  of  an  accurately  Gaussian  component,  however,  requires  a  good  deal  of  data  manip¬ 
ulation  and  sets  aside  another  component  of  complex-signal  scintillation. 

Figure  5  also  illustrates  (with,  in  fact,  unusual  fidelity)  the  consistent  result  that 
phase  scintillations  are  well  described  by  a  normal  distribution.  It  demonstrates,  further, 
that  the  two-component  model  retains  this  essential  behavior,  yielding  a  phase  pdf  only 
slightly  less  efficacious  than  the  log-normal  model  and  far  better  than  the  generalized- 
Gaussian  hypothesis.  In  contrast,  the  generalized-Gaussian  model  provides  no  mechanism 
for  describing  phase  beyond  a  modulo-2ir  definition. 

A  strictly  formal  test  of  the  phase  pdf  arising  from  the  generalized-Gaussian  hypoth¬ 
esis  would  be  to  compare  the  pdf  with  a  histogram  of  phase  values  for  which  2ir  ambiguities 
have  not  been  resolved.  To  do  so  in  the  present  work,  however,  would  have  been  to  treat 
the  data  differently  for  different  hypotheses;  it  was  decided,  instead,  to  treat  the  data 
identically  in  every  case  and  to  accept  poor  phase  fits  as  a  quantitative  statement  of  the 
qualitatively  obvious  limitation  of  generalized-Gaussian  statistics  in  the  presence  of  large 
phase  excursions. 

In  performing  the  tests,  it  was  supposed  that  the  phase  pdf  repeated  itself  in  every 
2tt  range  for  which  phase  values  were  observed.  Normalization  was  carried  out  to  match  the 
area  under  the  pdf  with  that  under  the  observed  histogram  over  the  range  of  bins  containing 
values.  An  alternative  would  have  been  to  match  the  area  under  the  pdf  in  the  range  -ir  to 
+tt  with  that  under  the  entire  histogram,  which  also  would  have  produced  large  x2  mismatches. 

A  further  point  demonstrated  in  Figure  5,  et  sequel,  is  the  efficacy  of  the  Nakagami-m 
distribution  in  describing  intensity  scintillation.  As  was  discussed  in  Section  III  A  1 
and  illustrated  in  Tables  3  and  4,  this  was  one  of  the  most  consistent  results  of  hypoth- 
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component  (top)  and  the  logarithm  of  the  focus  component  (bottom)  of  the 
L-band  signal  shown  at  the  bottom  of  Figure  7,  compared  with  normal  pdf's 
computed  from  the  measured  means  and  variances. 
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esis  testing.  Figure  5  points  up  particularly  well  that  a  bivariate  distribution  re¬ 
ducing  to  the  product  of  the  Nakagami-m  distribution  for  intensity  and  the  normal  distrib¬ 
ution  for  phase  in  cases  of  no  correlation  between  the  variates,  would  provide  a  suitable 
model  for  complex-signal  scintillation. 

Figure  7,  however,  illustrates  the  fact  (particularly  at  UHF)  that  such  correlation 
is  a  salient  feature  of  scintillating  signals.  It  may  well  be  unimportant  in  many  appli¬ 
cations,  but  it  cannot  be  ignored  in  a  general  statistical  description  of  scintillation. 
Furthermore,  it  can  be  calculated  as  a  function  of  Fresnel  distance  fFremouw  and  Rino, 

1978).  That  the  two-component  hypothesis,  which  does  account  for  correlation,  provides 
an  efficacious  model,  also  is  illustrated  in  Figure  5,  et  sequel.  The  two-component  model 
also  provides  some  intuitive  basis  for  understanding  signal-statistical  behavior,  although 
at  the  expense  of  requiring  six  parameters  (the  "six  sigmas")  for  its  specification. 

We  have  dealt  only  with  raw  \2  values  for  assessing  efficacy  of  the  four  models.  In 
general,  the  raw  \2  values  are  consistent  with  the  impressions  one  gains  from  inspection  of 
graphical  histogram  fits.  For  hypothesis  testing  in  the  strict  sense,  one  would  want  to 
establish  levels  of  confidence,  which  can  be  done  for  individual  data  sets  but  not  for  the 
mean  \2  values  appearing  in  Table  3. 

Even  for  individual  data  sets,  subtle  questions  arise  regarding  the  number  of  degrees 
of  freedom  lost  in  computing  the  intensity  and  phase  pdf's  on  the  basis  of  complex-signal 
moments  (for  the  generalized  Gaussian  and,  particularly,  the  two-component  hypotheses). 
Suffice  it  to  say  that  a  wide  range  of  confidence  levels  (from  0.01%  to  99.5%)  has  been 
calculated  on  the  basis  of  reasonable  assumptions  about  degrees  of  freedom  lost  and  that 
the  rankings  of  the  various  hypotheses  have  remained  the  same  as  those  established  by  means 
of  the  raw  x2  values. 

B.  GEOMETRICAL  CONTROL  OF  THE  RATIO,  S,  /ci 

4  tp 

l.  Relevant  Theoretical  Expressions 

In  Section  III  A,  we  explored  the  variation  of  complex-signal  statistics  with  observing 
frequency,  scattering  strength,  and  receiver  location  but  postponed  a  detailed  description 
of  geometrical  factors.  The  purpose  of  this  section  is  to  provide  such  a  description,  not 
in  terms  of  the  statistical  distributions  but  in  terms  of  the  two  most  relevant  moments 
that  control  them:  namely,  the  phase  and  intensity  scintillation  indices,  o  and  S,  (=ra  ~). 
One  of  the  first  features  of  complex-signal  statistics  noted  in  Wideband  data  was  a  size¬ 
able  difference  in  the  ratio  of  S,  to  o,  when  measured  at  Stanford,  Poker  Flat,  and  Ancon 
(Fremouw  and  Livingston,  1976).  Later  it  was  found  that  both  equatorial  stations,  Ancon 
and  Kwajalein,  consistently  show  a  higher  S.  /a  ratio  than  does  Poker  Flat,  for  nonsaturated 
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intensity-scintillation  conditions  (Fremouw,  1978a). 
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Rino  has  developed  weak-scatter  expressions,  based  on  the  phase-screen  propagation 
model,  that  are  suitable  for  interpreting  both  phase  (Rino  and  Fremouw,  1977)  and  inten¬ 
sity  (Rino  and  Matthews,  1978)  scintillation  records  in  terms  of  a  fairly  general  descrip¬ 
tion  of  irregular  plasma  structure.  We  employ  here  his  expressions  for  a  and  ,  finding 

that  the  observed  differences  in  their  ratio  may  be  explained  in  terms  of  scattering- 

geometry  control  in  the  presence  of  an  outer  scale  that  falls  outside  our  data  window.  We 

suggest  that  the  same  effect  causes  the  latitudinal  and  Fresnel-zone-size  dependences  of 
best-iit,  signal-statistical  models  noted  in  Table  3. 

Recalling  that  the  detrending  procedure  employed  in  Wideband  data  processing  behaves 
like  a  sharp-cutoff,  high-pass  filter  (Fremouw,  1978b),  we  may  write  the  square  of  the 
phase  scintillation  index  succinctly  as 

2T  f(,_2v) 

-  (26) 


where  T  =  power-spectral  density  of  phase  at  1  Hz, 
f  “  detrender  cutoff  frequency, 

and  v  =  power-law  spectral  index  describing  the  spatial  spectrum  of  scattering 
irregularities  (obtained  with  resonable  accuracy  as  half  the  measured 
index  of  the  phase  power  spectrum). 

Rino's  formulation  relates  the  measured  strength  of  the  phase  spectrum  to  ionospheric 
and  geometrical  parameters  as 


T  =  (r2A2L  sec9)  G  Cg 


(2tt  )  2  v+1  p  (v+£) 


where  r  =  classical  electron  radius, 

e 

\  =  radio  wavelength, 
and  6  =  incidence  angle. 


The  gamma  functions  arise  from  normalizing  the  three-dimensional  ionospheric  spectrum  to 
^(AN)2^>,  and  the  main  geometrical  factor  is 

G  =  .  -  (28 

VAC  -  B2/4  ros'l 


where  a  and  b  are  irregularity  axial  ratios  along  and  across  the  geomagnetic  field,  respec¬ 
tively,  and  A,  B,  and  C  are  geometrical  parameters  given  in  Eq ,  (41)  of  Rino  and  Fremouw 
(1977).  The  strength  of  the  scattering  structure  is  described  as 
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8tt 


3/2 


r  (v~H)  2  v— 2 

r(v-i)  q0 


<(an2)> 


(29) 


which  is  (to  within  a  factor  ab,  which  Rino  chose  to  incorporate  in  G)  numerically  equal  to 
the  strength  of  the  "isotropized"  three-dimensional  spectrum  at  a  wave  number  of  1  radian/m. 
The  spectrum  is  supposed  to  be  cut  off  at  a  very  large  outer  scale,  q 

The  remaining  quantity  in  Eq.  (27)  is  the  "effective  velocity",  Ve ,  which  is  the  speed 
with  which  the  radio  line  of  sight  cuts  across  contours  of  equal  correlation  in  the  iono¬ 
sphere  (or  alternatively  with  which  contours  of  phase  correlation  scan  across  the  re¬ 
ceiver).  The  effective  velocity  accounts  for  the  effect  that  anisotropy  of  the  irregular¬ 
ities  has  on  the  conversion  of  spatial  structure  to  temporal  fluctuations;  it  is  related 
to  scan  velocity  in  isotropic  coordinates  by  means  of  Eq.  (12)  of  Rino  and  Matthews  (1978). 

The  square  of  the  intensity  scintillation  index  developed  in  Rino  and  Matthews  (1978) 
is 


S2  =  4  (r2X2  L  sec  9)  F  C 
4  e  s 


-r  (■*=?*)  cos  [u(l-2v)/4] 


2u  2 


(5-2 v) / 2 


in  which  diffraction  is  described  primarily  by  the  Fresnel-zone  parameter 
AzsecO 


Z  = 


4  IT 


(30) 


(31) 


where  z  is  the  effective  "reduced  height"  (including  correction  for  wavefront  curvature)  of 
the  irregularities.  The  major  theoretical  result  presented  in  Rino  and  Matthews  (1978)  is 
evaluation  of  the  geometry-dependent  Fresnel  filter  factor,  F,  as 


where  A"  and  C"  are  geometrical  parameters  derived  from  A,  B,  and  C  by  means  of  a  coordinate 
rotation  and  where  ?Fj  Is  the  Gaussian  hypergeometric  function.  Rino  and  Matthews  (1978) 
combined  Eqs.  (27)  and  (30)  into  the  following  succinct  description  of  the  dependence  of  the 
intensity  scintillation  index  on  the  strength  of  the  phase  perturbation  produced  by  the 
ionosphere : 


=  4C(v) _ - _ -  - 

’4  y  (Zv-1)  G 


(33) 


The  dependence  of  intensity  scintillation  on  the  shape  of  the  ionospheric-irregularity 
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power  spectrum  is  a  bit  complicated,  witn  i  In-  spectral  index  appearing  in  Eq .  (33)  in  four 
different  factors,  including 


-(27i)2vr(i-v)r(A+v)  cosT"  '  1  1 

C  (v)  = - - - — i J 


\  tt  r (u)  2 


(-2V)/? 


(34) 


-2 

For  a  one-dimensional  irregularity  spectrum  of  k  ,  as  is  often  observed  in  situ  (Dyson, 
McClure,  and  Hanson,  1974),  v  =  1.5.  Unfortunately,  Eq.  (34)  takes  on  the  indeterminate 
value  0/0  in  this  case  [as  does  the  factor  in  curly  brackets  in  Eq.  (30)]. 

However,  by  applying  recurrence  and  reflection  formulas  for  the  gamma  function 
(Abramowitz  and  Stegun,  1964)  and  a  trigonometric  identity,  one  may  show  that  Eq.  (34)  is 
equivalent  to 


C(v)  = 


\  n  (2tt)2'~  (cos  cos  ~  +  sin  ~  sin  -^p) 
2(5  2v^2p(v)cos  (ttv) 


(35) 


which,  by  means  of  some  additional  trigonometric  relations,  may  be  simplified  to 
23vtt2v  cos  ("ft  (^  -  i)l 

C(v)  =  - —  - L - J  (36) 

— 4v  2 /rr  T  (v) cos  (vtv) 

A  single  application  of  L'Hopital's  rule  then  shows  that 


Lim  C(v)  =  4tt3  (37) 

v->-l .  5 

Figure  10  is  a  plot  of  Eq.  (36)  and  (37),  showing  that  C(v)  mitigates  somewhat  the  tendency 

for  to  increase  with  decreasing  spectral  index  via  the  Fresnel  filter  factor,  F,  and 

the  factor  V  1  2v  in  Eq.  (33),  as  doe^  the  factor  7.V 
e 

2.  Application  to  Investigation  of  the  Scintillation-Index  Ratio 

Eq.  (33)  gives  the  intensity-scintillation  index  in  terms  of  the  phase  power-spectral 
density,  T.  To  obtain  the  intensity  variance  per  unit  phase  variance  experienced  by  a  sys¬ 
tem  with  a  high-pass  cutoff,  f  ,  we  divide  Eq.  (33)  by  Eq.  (26)  to  obtain 


[2  (2v-l)C(v) 


(38) 
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Various  observers  have  reported  differences  in  the  spectral  index,  v,  inferred  from  scin¬ 
tillation  measurements  made  at  different  locations  (Whitney  and  Basu,  1977;  Rino  et  al, 
1978;  Rino  and  Matthews,  1978),  which  could  produce  differences  in  the  scintillation-index 
ratio.  To  the  extent  that  the  spectral  index  displays  a  systematic  latitude  dependence, 
however,  it  appears  to  be  in  the  sense  of  smaller  values  of  v  in  auroral  than  in  equatorial 


regions.  This  is  opposite  the  sense  required  to  explain  the  observed  S.  /a, 
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study  here. 


behavior  under 


The  other  obvious  candidate  for  explaining  the  observed  behavior  is  the  possibility 


of  different  effective  heights  for  the  relevant  striations  at  different  latitudes,  which 


would  affect  Eq.  (38)  directly  through  Z.  Rino  and  Matthews  (1978)  have  pointed  out. 


however,  that  this  dependence  is  complicated  by  the  indirect  effect  of  irregularity  height 
on  Ve  (by  way  of  the  scan  velocity  for  a  given  satellite). 

In  the  work  reported  in  this  section,  we  held  v  and  the  irregularity  height  constant 
in  order  to  investigate  the  role  that  factors  having  to  do  with  striation  shape  and  ex¬ 
periment  geometry  might  play  in  establishing  the  Sj,  /o^  ratio.  We  then  performed  comparative 
calculations  for  VHF  data  sets  from  two  Wideband  passes  -  Poker  Flat  Pass  6-39,  observed 
on  28  July  1976,  and  Kwajalein  Pass  22-13,  observed  on  28  October  1977.  Figure  11  illus¬ 
trates  the  reason  for  selecting  these  two  samples,  which  contain  60  seconds  of  data  from 
near  1020  GMT  (0020  local  time)  and  85  seconds  of  data  from  near  1228  GMT  (0028  local  time), 


respectively. 

The  figure  shows  plots  of  versus  for  VHF  data  sets  from  all  four  Wideband  sta¬ 
tions  for  which  dual-detrend  analysis  has  been  carried  out.  The  data  sets  were  carefully 
selected  to  include  a  wide  range  of  observing  geometries  and  scintillation  levels  and  to 
provide  reasonable  statistical  stationarity.  The  two  data  sets  marked  on  the  figure  were 

chosen  for  this  investigation  of  S.  la.  because  they  are  representative  of  results  at  an 

4  <p 

equatorial  and  an  auroral-zone  station  under  conditions  of  essentially  identical,  moderate 
intensity  scintillation. 

To  investigate  the  clear-cut  difference  between  behavior  at  Poker  Flat  and  at  the  two 
equatorial  stations,  the  known  observing  geometry  associated  with  the  two  data  points  iden¬ 
tified  in  Figure  11  was  employed,  along  with  several  irregularity  models,  to  calculate 
various  factors  appearing  in  Eq.  (38).  Finally,  the  factors  were  combined  to  calculate  the 
scintillation-index  ratio  itself,  the  results  of  which  are  shown  in  Figure  12. 

The  first  model  employed  was  that  of  axially  symmetric  striations  aligned  along  the 
geomagnetic  field.  For  this  model,  the  single  axial  ratio  was  stepped  from  unity  (isotropic 
irregularities)  to  40:1,  as  shown  on  the  steep  solid  curve  in  the  figure.  Values  of  S 
were  calculated  from  the  model  and  from  the  known  ionospheric  incidence  angle  and  geomagnetic 
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POKER  FLAT  PASS  6-39 

The  ratio  of  intensity  scintillation  index,  S,,  to  phase  scintillation  index, 
oa,  for  a  50-second  data  set  from  Poker  Flat  (horizontal  axis)  and  an  85- 
second  data  set  from  Kwajalein  (vertical  axis).  Isolated  point  at  (0.24, 
1.10)  represents  observed  data.  Other  points  are  calculated  from  various  rod 
and  sheet  irregularity  models.  Numbers  by  points  are  axial  ratios.  (All 
sheet  models  have  two-dimensional  isotropy;  that  is,  10  means  10:10:1  and  40 
means  40:40:1.) 
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heading  of  the  propagation  vector  and  scan  velocity  of  the  line  of  sight,  for  each  of  the 
two  data  sets.  The  results  for  the  two  samples  were  then  plotted  on  a  common  grid,  as 
were  the  observed  scintillation-index  ratios  (isolated  point  at  0.24,  1.10). 

Next,  sheet-like  irregularity  models  were  invoked  and  the  calculations  repeated. 

Three  such  models  were  employed:  sheets  layered  like  onion  skins  along  L  shells,  such  as 
have  been  observed  from  Poker  Flat  (Fremouw  et  al,  1977)  and  possibly  at  a  middle  latitude 
in  the  southern  hemisphere  (Hajkowicz,  1978);  vertical  sheets  lying  in  the  magnetic  mer¬ 
idian  plane;  and  sheets  oriented  at  45°  to  both  the  meridian  plane  and  the  L  shells,  tilted 
from  upper  west  to  lower  east  as  some  backscatter  plumes  appear  to  be  (Fremouw  and 
Lansinger,  1977;  Tsunoda,  Baron,  and  Owen,  1978).  For  each  sheet  model,  the  calculation 
was  performed  for  axial  ratios  of  10:10:1  and  40:40:1. 

The  first  point  to  be  noted  from  Figure  12  is  that  all  models  but  one  (40:40:1  sheets 

aligned  on  L  shells)  produced  a  greater  S, /o,  ratio  for  Kwajalein  than  for  Poker  Flat. 

4  <p 

Thus,  one  is  led  immediately  to  suspect  that  something  inherent  in  the  geometry  produces 
the  observed  difference.  It  happens  that  the  ionospheric  incidence  angle,  6,  is  larger 
for  the  Kwajalein  data  set  (54.5°)  than  for  the  set  from  Poker  Flat  (17.9°).  This  fact 
alone  is  sufficient  to  produce  a  slightly  greater  scintillation-index  ratio  in  the  K22-13 
data  set  than  in  the  PF6-39  set,  as  seen  in  the  point  calculated  for  isotropic  irregular¬ 
ities  (1:1  rods).  It  is  not  sufficient,  however,  to  reproduce  the  observed  data  point. 
Indeed,  only  five  of  the  12  models  employed  produced  calculated  points  lying  on  or  above  a 
line  between  the  origin  of  Figure  11  and  the  observed  data  point. 

There  is  no  reason  to  suppose  a  pAXOAx  that  the  same  irregularity  model  is  appropriate 
at  both  Kwajalein  and  Poker  Flat.  As  an  aside,  it  is  interesting  to  note  that  few  of  the 
models  invoked  produce  enough  diffraction  to  yield  S^/o^  ratios  as  large  as  those  observed 
at  either  station  (for  the  assumed  values  of  v  and  irregularity  height),  although  several 
come  close  at  Poker  Flat.  Highly  elongated  rod-like  striations,  however,  have  no  difficulty 
producing  the  large  /o^  ratio  observed  at  Kwajalein. 

The  main  thrust  of  this  analysis  is  to  attempt  to  isolate  the  factor  that  dictates  the 
large  disparity  in  scintillation-index  ratio  observed  at  the  two  stations.  The  final  two 
figures  present  a  comparison  between  two  relevant  factors.  Figure  13  shows  calculated  re¬ 
sults  for  the  ratio  of  the  Fresnel-dif fraction  filter  factor,  F,  to  the  phase-scintillation 
geometric  enhancement  factor,  G.  This  ratio  contains  most  of  the  "static"  effect  of  dif¬ 
fraction.  That  is,  it  is  independent  of  source  motion  and  therefore  of  temporal  filtering 
operations  such  as  detrending. 

The  key  point  to  note  in  Figure  13  is  that  all  models  produce  calculated  points  lying 
on  (isotropic  only)  or  below  (all  other  models)  a  45°  radial  from  the  origin  of  the  calcu- 
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Figure  13.  Calculated  values  of  the  main  geometrical  factor  controlling  "static"  dif¬ 
fraction:  i.e.,  the  ratio  of  the  Fresnel-dif f ract ion  filter  factor,  F,  to 
the  phase-scintillation  enhancement  factor,  G. 


lation  grid.  That  is,  all  of  the  models  invoked  produced  more  diffraction  at  Poker  Flat 
than  at  Kwajalein  (aside  from  sec  9  factors  that  are  not  included  in  F  and  G) .  While 
not  anticipated,  this  fact  probably  results  from  the  element  of  alignment  along  the  geo¬ 
magnetic  field  (which  is  horizontal  at  Kwajalein  and  highly  inclined  at  Poker  Flat)  that 
is  inherent  in  all  the  models.  The  main  point  is  that  effects  of  static  geometry  are  in 
the  opposite  sense  from  that  needed  to  explain  the  observed  disparity  in  scintillation- 
index  ratio. 

The  situation  reverses  abruptly  when  one  considers  source  motion  and  the  effect  of 
detrending  on  phase  scintillation,  as  demonstrated  in  Figure  14.  The  ratio  f  /V  is  the 
spatial-filter  cutoff  (in  units  of  m  imposed  by  the  detrender  on  the  isotropically  ex¬ 
pressed  nonisotropic  spatial  spectrum  of  phase  fluctuations.  In  effect,  it  translates  the 
detrend  filter  not  only  from  time  to  space  but  also  from  the  actual  scan  direction  to  an 
euqivalent  scan  across  the  narrow  dimension  of  the  scattering  irregularities  (i.e.,  the 
dimension  having  size  unity  in  an  a:b:l  model  with  a  and  b  1  1).  Figure  14  contains  plots 
of  (fQ /Ve ) 2v-l  wlth  v  set  to  1.5.  Thus,  10  8  corresponds  to  a  spatial-filter  cutoff  of 
10  km  and  10  6  to  a  cutoff  of  1  km. 

What  we  see  in  Figure  14  is  that  the  detrender  cutoff  is  at  higher  spatial  frequencies 
in  the  isotropized  spectrum,  for  all  models,  at  Kwajalein  than  at  Poker  Flat.  (That  is, 
all  calculated  points  lie  above  a  45°  radial.)  The  small  difference  for  isotropic  (1:1) 
irregularities  results  from  simple  projection  effects  having  to  do  with  the  different  look 
angles  to  the  satellite  for  the  two  data  sets. 

The  dotted  diagonal  line  in  the  lower  portion  of  the  figure  corresponds  to  the  isolated 
observational  data  point  in  Figure  12.  That  is,  a  model  that  produces  a  point  on  that  line 
satisfies  both  observations.  Again,  there  is  no  reason  to  expect  a  single  model  to  satisfy 
both  observations;  the  point  is  that  all  models  except  nearly  isotropic  ones  result  in 
effective  spatial-filter  cutoffs  due  to  detrending  that  can  readily  explain  the  large  dis¬ 
parity  in  the  scintillation-index  ratio  for  the  two  data  sets. 

3.  Discussion 

It  is  intuitively  understandable  that  detrender  action  can  affect  the  scintillation- 
index  ratio  because  it  directly  sets  the  high-pass  cutoff  of  phase  scintillation,  whereas 
the  diffraction  process  itself  performs  that  function  for  intensity  scintillation,  by  means 
of  Fresnel  filtering.  Thus,  it  was  suspected  a  priori  that  the  detrender,  coupled  with 
irregularity  anisotropy  and  line-of-sight  scanning,  might  play  a  role  in  the  differences  ob¬ 
served  in  the  scintillation-index  ratio  at  different  stations.  Figure  14  shows  that  the 
suspected  effect  is  in  the  correct  sense  to  explain  the  observation  and  is  very  strong. 

Indeed,  virtually  all  field-aligned  models  produce  sufficiently  anisotropic  effective 
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POKER  FLAT  PASS  6-39 

Figure  14.  The  factor  (f0/Ve)2v  1  which  arises  from  spatial  filtering  of  anisotropic  ir¬ 
regularities  by  means  of  the  Wideband  detrender.  Dashed  diagonal  line  in 
lower  portion  of  figure  corresponds  to  observed  data  point  in  Figure  12. 
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spatial-filter  cutoffs  that  other  factors  must  be  operating  to  prevent  even  greater  dis¬ 
parity  in  the  index  ratio  observed  at  equatorial  stations  and  at  an  auroral-zone  location 
(for  a  high-inclination  satellite  such  as  Wideband).  One  such  factor  is  the  static  dif¬ 
fraction  factor  represented  by  F/G.  Another  may  be  less  steep  irregularity  spectra  at 
auroral  than  at  equatorial  locations. 

The  result  of  this  analysis  underscores  the  central  role  played  by  a  steep,  effect¬ 
ively  low-pass  spatial  spectrum  of  irregularities  in  dictating  the  signal  statistics  of 
scintillation.  That  is,  since  the  outer  scale  of  ionospheric  structure  lies  outside  the 
spectral  window  of  structures  that  produce  scintillation,  an  effective  outer  scale  is  im¬ 
posed  by  the  act  of  phase  detrending. 

The  pertinent  point  is  that  the  outer  scale  imposed  is  geometry-dependent.  The  result 
is  that  the  Fresnel  distance  (ratio  of  Frensel-zone  size  to  effective  outer  scale)  depends 
strongly  on  the  direction  of  scan  relative  to  elongation  axes  of  refractive-index  structure. 
Detrending,  therefore,  reintroduces  the  concepts  of  "near"  and  "far"  zones  which  appear  to 
be  lost  in  a  strict  power-law  spectral  environment.  The  relative  efficacy  of  log-normal 
(near-zone)  and  generalized-Gaussian  (far-zone)  complex-signal  statistics  at  various  sta¬ 
tions  and  different  frequencies  may  be  understood  in  these  terms. 

In  closing  this  section,  we  stress  that  detrending  effects  are  not  to  be  viewed  merely 
as  artifacts  of  the  Wideband  experiment.  Such  effects  arise  in  all  experiments  and  oper¬ 
ational  systems,  and  they  are  important  factors  in  controlling  the  quantitative  relationship 
between  phase  and  intensity  scintillation  and,  therefore,  the  statistics  of  complex-signal 
scintillation. 
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IV.  CONCLUSION 


The  Wideband  experiment  offers  a  good  deal  of  insight  into  the  behavior  of  signals 
chat  have  undergone  narrow-angle  forward  scattering  in  a  random  medium.  From  the  present 
work,  one  may  state  several  conclusions  about  the  underlying  statistics  by  which  such 
signals  may  be  described.  First,  both  of  the  classical  views  —  ie,  "multiplicative"  and 
"additive"  (Strohbehn,  Wang,  and  Speck,  1975)  —  contain  elements  of  applicability  to 
complex-signal  statistics. 

The  multiplicative  view,  however,  is  appropriate  only  for  describing  phase  fluctuations 
and  weak  intensity  scintillation.  Phase  is  nearly  always  normally  distributed  without 
regard  to  the  strength  of  scatter  or  size  of  the  Fresnel  zone.  Thus,  its  first-order  sta¬ 
tistics  are  altered  very  little  by  post-scattering  propagation,  caveats  being  that 
sufficient  signal  strength  be  maintained  that  phase  remains  defined  and  that  the  signal  be 
sampled  frequently  enough  that  phase  ambiguities  are  resolved. 

Regarding  intensity,  the  present  analysis  consistently  supports  the  additive  view. 

That  is,  the  log-normal  distribution  for  intensity  is  efficacious  only  in  circumstances 
where  there  is  little  difference  between  it  and  other  distributions  (small  Fresnel  zone 
and/or  weak  scatter)  or,  occasionally,  in  situations  where  statistical  significance  is 
limited  by  lack  of  stationarity. 

A  strong  conclusion  can  be  made  in  favor  of  the  Nakagami-m  distribution  for  describing 
intensity  scintillation.  Often  regarded  as  equivalent  to  the  Rice  distribution  over  the 
range  of  m  usually  encountered  under  scintillation  conditions,  the  Nakagami  distribution 
actually  is  an  approximation  to  a  more  general  intensity  distribution  consistent  with  the 
generalized  Gaussian  (additive)  hypothesis.  (cf.  Section  3.1  (b)  of  Nakagami,  I960.)  The 
Rice  distribution  was  not  tested  specifically  in  this  work,  but  the  Nakagami  distribution 
was  found  to  be  consistently  efficacious  over  the  range  of  m  in  which  the  two  are  somewhat 
different . 

The  conclusions  quoted  above  lead  to  a  suggestion  that  complex-signal  statistics  could 
be  described  very  well  by  a  bivariate  distribution  for  intensity  and  phase  that  would  re¬ 
duce  to  the  product  of  the  Nakagami-m  distribution  and  the  normal  distribution,  respectively, 
when  correlation  between  the  variates  is  zero.  In  general,  however,  non-zero  correlation 
is  observed  experimentally  and  is  calculable  from  propagation  theory.  It  is  suggested  that 
an  inductive  procedure-  starting  from  the  Nakagami  and  normal  distributions  and  invoking  the 
necessary  bivariate  normalization  might  lead  to  a  simple  and  efficacious  model  for  complex- 
signal  statistics  that  would  adequately  describe  correlation  between  intensity  and  phase. 

At  present,  only  two-component  hybrid  models  are  able  to  provide  a  full  description  of 
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complex-signal  statistics.  The  two-component  model  of  Fremouw,  Rino,  and  Livingston 
(1976),  tested  in  this  work,  provides  a  valid  description  of  complex-signal  behavior  for 
all  relevant  regimes  of  scattering  strength  and  Fresnel-zone  size.  It  does  so,  however, 
at  the  expense  of  requiring  six  parameters  for  its  specification.  A  model  based  on  the 
Nakagami  and  normal  distributions,  as  suggested  in  the  foregoing  paragraph,  would  require 
only  three  parameters.  For  applications  not  requiring  the  joint  statistics  of  (correlated) 

intensity  and  phase,  the  Nakagami-m  and  normal  distributions  can  provide  very  efficacious 

_2 

descriptions  on  the  basis  of  a  single  measured  parameter  each  (m  =  S.  and  o  ,  respectively) 

4  <P 

The  two-component  model  yields  insight  into  the  physical  processes  governing  signal 
behavior,  in  that  it  provides  a  basis  for  isolating  effects  adequately  described  by  means 
of  ray  optics  and  effects  dominated  by  diffractive  scatter.  In  the  model,  these  components 
are  taken  to  obey  log-normal  and  generalized  Gaussian  statistics,  respectively.  It  is  only 
in  this  context  that  generalized  Gaussian  statistics  provide  a  (partial)  description  of 
complex-signal  scintillation.  Thus,  while  intensity  scintillations  are  well  described  by 
distributions  based  on  the  generalized-Gaussian  hypothesis,  as  was  found  by  Rino,  Livingston 
and  Whitney  (1976)  and  as  is  totally  consistent  with  the  efficacious  performance  of  the 
Nakagami-m  distribution  found  in  the  present  work,  that  hypothesis  is  not  supported  by, 
nor  does  it  adequately  describe,  complex-signal  scintillation. 

In  a  strict  sense,  the  data  employed  in  this  analysis  describe  signal  behavior  only  in 
the  "intermediate"  zone,  in  which  the  Fresnel-zone  radius  lies  between  the  inner  and  outer 
scales  of  a  medium  described  by  a  power-law  spatial  spectrum.  In  principal,  the  inner 
scale  was  effectively  zero  (smaller  than  the  spatial  sampling  size  and/or  lying  at  a  point 
in  the  spatial  spectrum  that  was  dominated  by  noise)  and  the  outer  scale  was  effectively 
infinite  (larger  than  the  largest  scale  admitted  in  a  data  segment).  Application  of  a  de¬ 
trending  filter  to  the  measured  phase,  however,  had  the  effect  of  imposing  an  outer  scale 
and  allowed  investigation  of  near-zone  and  far-zone  behavior.  Such  investigation  was  facil¬ 
itated,  also,  by  the  multi-frequency  nature  of  the  Wideband  experiment. 

Temporal  detrending  in  the  presence  of  highly  anisotropic  scatterers  results  in  dif¬ 
ferent  imposed  outer  scales  for  different  observing  geometries  —  in  particular,  for  those 
encountered  in  observation  of  a  high-inclination  satellite  from  stations  located  at  low, 
middle,  and  high  latitudes.  Such  detrending,  which  is  imposed  in  some  fashion  in  nearly  all 
measurements  and  applications,  has  a  strong  effect  on  the  observed  ratio  of  S,  to  o.  and, 
thereby,  on  the  statistical  behavior  of  the  complex  signal.  This  strong  geometrical  control 
over  the  relationship  between  intensity  and  phase  statistics  must  be  taken  into  account  in 
any  application  dealing  with  total  variances  of  refractive-index  structure  and/or  signal 
fluctuations  (as  contrasted  with  power  spectral  densities  of  these  same  random  variables). 
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